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ABSTRACT 

We present the results of an extensive exploration of the five-dimensional parameter 
space of the minimal SU(5) supergravity model, including the constraints of a long enough 
proton lifetime (r p > 1 x 10 32 y) and a small enough neutralino cosmological relic density 
(Q x /iq < 1). We find that the combined effect of these two constraints is quite severe, 
although still leaving a small region of parameter space with rrig^q < 1 TeV. The allowed 
values of the proton lifetime extend up to r p ~ 1 x 10 33 y and should be fully explored 
by the Super Kamiokande experiment. The proton lifetime cut also entails the following 
mass correlations and bounds: mn ^ 100 GeV, m x « \m x ^ ~ 0.15mg, m x o « m x +> 
and m Y < 85 (115) GeV, m v o v + < 165 (225) GeV for a 3 = 0.113(0.120). Finally, the 
combined proton decay and cosmology constraints predict that if rrih ^ 75 (80) GeV then 
m + ^ 90 (110) GeV for as = 0.113 (0.120). Thus, if this model is correct, at least one of 
these particles will likely be observed at LEPII. 
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1. Introduction 

The minimal SU(5) supergravity model [Q has recently passed the simplest possible 
consistency check, namely the unification of the gauge couplings at an energy scale Mjj ~ 
10 16 GeV [p]. This check depends only in sub-leading order on the masses of the light 
and heavy particles in the theory, and as such provides weak constraints on the various 
model parameters || . On the other hand, the requirement that the dimension- five-induced 
proton decay operators be within current experimental bounds provide rather stringent 
constraints on all sectors of the theory [p|,|6|,[7i|8|,|9|,|T0| . Recently two of us ( JL and DVN) with 



A. Zichichi |TT| , studied a representative set of points in parameter space which satify the 
proton decay bound and applied to these the cosmo logical requirement O x /iq < 1, where Vt x 
is the relic abundance of the lightest neutralino (also the lightest supersymmetric particle 
(LSP) and which is assumed to be stable) and 0.5 < ho < 1 is the Hubble parameter. 
We found that the cosmological constraint was grossly violated for these points. It was 
also noted that there may still exist cosmologically allowed regions for sufficiently small 
values of mt. In this paper we present a systematic exploration of the five-dimensional 
parameter space of the model, which corroborates the previous indicative results. We 
identify the small regions of parameter space which satisfy both constraints and show 
that the cosmologically acceptable region should be almost fully testable at LEPII and/or 
SuperKamiokande. 

The minimal SU(5) supergravity model consists of the Standard Model fields (with two 
Higgs doublets) and their superpartners, plus the heavy GUT fields in the form of colored 
Higgs triplet fields (if, H) and gauge and Higgs bosons in the 24 representation of SU (5). 
Universal soft-super symmetry breaking at the scale Mjj is described by three parameters: 
m 1( / 2 ,mo, A. The SU(5) gauge symmetry entails gauge coupling (a± = «2 = «3 = «5 
at Mjj) and Yukawa coupling (At, = A T at Mu) unification. In our calculation all these 
parameters are run from Mjj down to low energies (Mz) as prescribed by the relevant 
renormalization group equations. Finally, at the scale Mz the full one-loop effective po- 
tential is minimized and the one-loop corrected Higgs boson masses are obtained. The 
final parameter list is: tan/3, m t , TO1/2, £o = m o/ m i/2, £,A = ^-/ m i/2 ? an d the sign of 
the Higgs mixing term (n). This set of parameters should be constrasted with the 21- 
dimensional parameter space of the usual MSSM. For recent analyses of this type, see 
Refs. [0,|T3|Mff5T|T6|. 
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2. Standard Constraints 



We impose the following set of 'standard constraints' on the parameter space of the 
model (see Ref. jnjj for a detailed discussion): 

(i) one-loop radiative electroweak symmetry breaking; 

(ii) perturbative unification (which implies mt ^ 190 GeV and tan/3 < 50); 

(iii) ml r > 0; 

<?>» 

(iv) a neutral and colorless LSP (i.e., v or x = Xi); 

(v) experimental bounds on rn- x +, rrif, m g , m^, m t , T inv ; 

(vi) m~" e , m g < 1 TeV, motivated by naturalness considerations, or demanding testability 
of the model at the SSC; 

(vii) Xb(Mu) = X T (Mu), which in practice we use to predict m^m;,) for a given set of mt, 
tan/3, and a 3 values. 

We have explored the following hypercube of the parameter space: \i > 0, fx < 0, 
tan/3 = 2 - 10(2), m t = 100 - 160(5), £ = - 10(1), £ A = -f ,0, +f , and m 1/2 = 
50 — 300 (6), where the numbers in parenthesis represent the size of the step taken in that 
particular direction. (Points outside these ranges have little (a posteriori) likelihood of 
being acceptable.) Of these 92, 235 x 2 = 184, 470 points, « 25% passed all the standard 
constraints. In what follows we present our results as a collection of scatter plots, where a 
given pair of observables is plotted for each allowed point in the parameter space. 



3. Proton Decay 

In unified supersymmetric theories only the dimension-five-mediated proton decay 
operators are constraining. In calculating the proton lifetime we consider the typically 
dominant decay modes p — > u^^ T K + and neglect all other possible modes. Schematically 
the lifetime is given by M 



M H sin 2/3 1 ' 



fl+ytK 



(3.1) 



Here Mh is the mass of the exchanged GUT Higgs triplet which on perturbative grounds 
is assumed to be bounded above by Mh < 3Mjj [|6|,PJI0[]; sin2/3 = 2tan/3/(l + tan 2 /?), 
thus t p 'likes' small tan/3 (we find that only tan/3 < 6 is allowed); 1 + y tK represents 
the calculable ratio of the third- to the second-generation contributions to the dressing 
one-loop diagrams. An unkown phase appears in this ratio and we always consider the 
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weakest possible case of destructive interference. Finally / represents the sparticle-mass- 
dependent dressing one-loop function which decreases asymptotically with large sparticle 
masses. 

In Fig. 1 (top row) we show a scatter plot of (r p ,mg). The various 'branches' corre- 
spond to fixed values of Note that for £ < 3, t p < r^ xp = 1 x 10 32 y (at 90% C.L. flT7|). 
Also, for a given value of £o 7 there is a corresponding allowed interval in m§. The lower end 
of this interval is determined by the fact that r p oc l// 2 , and / ~ m x + ^ V m g( c + £o)' 
in the proton-decay-favored limit of \i Mjy; thus m g (c + £g) > constant. The upper 
end of the interval follows from our requirement m<j(oc rrig^c + £q) < 1 TeV. Statistically 
speaking, the proton decay cut is quite severe, allowing only about ~ 1/10 of the points 
which passed all the standard constraints, independently of the sign of fi. 

Note that if we take Mh = Mjj (instead of Mh = 3Mu), then r p — > ^r p and all 
points in Fig. 1 would become excluded. To obtain a rigorous lower bound on Mh, we 
would need to explore the lowest possible allowed values of tan/3 (in Fig. 1, tan/3 > 2). 
Roughly, since the dominant tan/3 dependence of r p is through the explicit sin 2/3 factor, 
the upper bound r p ^ 8x 10 32 y for tan (3 = 2, would become r p ^ lx 10 33 y for tan (3 = 1. 
Therefore, the current experimental lower bound on t p would imply Mh ^ Mjj- Note 
also that SuperKamiokande {r p xp ~ 2 x 10 33 y) should be able to probe the whole allowed 
range of r p values. 

The actual value of a^{Mz) used in the calculations («3 = 0.120 in all figures) has 
a non-negligible effect of some of the final results, mostly due to its effect on the value 
of Mh = 3M[/: larger values of as increase Mjj and therefore r p , and thus open up 
the parameter space, and viceversa. For example, for a.3 = 0.113(0.120) we get rrig ^ 
550 (800) GeV, £ > 5 (3), and r p < 4 (8) x 10 32 y. 



4. Cosmology 

We assume that the lightest neutralino is a stable particle, as expected in the minimal, 
.R-parity conserving model. The current cosmological observations of fio ^ 1 [0 and/or 
the inflation prediction f2o = 1 ||19|1 , lead us to impose the constraint fi x /io — 1- ^ n 
Fig. 2 we show the calculated valuesB of O x /iq versus m x , which show a noticeable dip at 



1 For a detailed discussion of the methods used to compute f2 x /ig see Ref. [pof . For computa- 
tions of f2 x /io in supergravity models with radiative electroweak symmetry breaking see [13,14,21]. 
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m x « \Mz due to s-channel Z-pole annihilation. Only ~ 1/6 of the points have O x /iq < 1. 
This result is not unexpected since proton decay is suppressed by heavy sparticle masses, 
whereas O x fog ^ s enhanced. Therefore, a delicate balance needs to be attained to satisfy 
both constraints simultaneously. The subset of cosmologically allowed points does not 
change the range of possible r p values (see Fig. 1 bottom row), although it depletes the 
const ant-£o ' branches ' . 

In Fig. 3 we show the calculated values of \x with and without the cosmological 
constraint. As expected, values of \x giving nearly pure gaugino x-compositions (above 
the \\x\= M2 = 0.3m g line) are cosmologically disfavored since in this limit the XX% an d 
XXh couplings are highly suppressed. Exceptions to this rule can be traced to m g values 
for which m x ~ \Mz, \ fnh when the annihilation cross section is enhanced by s-channel 
poles. Also, since grows with mt, larger values of mt tend to be disfavored as well. The 
results are quite similar for the 0:3 = 0.113 case. 

5. Particle Mass Correlations 

Since Fig. 3 shows that proton decay generally requires \i 3> Mw (and to a somewhat 
lesser extent also \x ^> M%), the lightest chargino will have mass ?^ x + ~ M% ~ 0.3m g , 
whereas the two lightest neutralinos will have masses m x kMi^ \M 2 and 
. Thus, within some approximation we expect 

m x ~ ^ m xa' m X2~ m x + ' rn x ^0.15m g . (5.1) 

The calculated values of these masses are shown in Fig. 4 (without imposing the cosmo- 
logical constraint). The approximate mass relations are quite accurate for ^ > 0, but more 
qualitative for [i < 0. Inclusion of the cosmological constraint basically just depletes the 
point density without affecting significantly the range of particle masses. The value of 
does not affect these mass relations either, although the particle mass ranges do change 

m y < 85 (115) GeV, m « + < 165 (225) GeV, for a 3 = 0.113 (0.120). (5.2) 

We also find that the one-loop corrected lightest Higgs boson mass (ra/J is bounded above 
by 

nih ^ 110 GeV, (5.3) 
4 



independently of the sign of //, the value of as, or the cosmological constraint. In Fig. 
5 we plot rrih versus m + which shows an experimentally interesting correlation if the 
cosmological constraints are imposed (bottom row). Indeed, 

m h > 75 (80) GeV m + < 90 (110) GeV, (5.4) 

for as = 0.113(0.120). Moreover, the correlations among the lightest chargino and neu- 
tralino masses in Eq. ( |5.1| ) imply analogous results for (to/i,to x o) and (m^,m x ), 



m h <; 80 GeV =>- m x o <> 90 (110) GeV, m x ^ 48 (60) GeV, (5.5) 

for «3 = 0.113 (0.120). These correlations can be understood in the following way: since 
we find that 3> Mz, then m/j ~ | cos 2f3\Mz + (rad. corr.). In the situation we consider 
here, we have determined that all of the allowed points for m g > 400 GeV correspond to 
tan/3 = 2. This implies that the tree-level contribution to rrih is ~ 55 GeV.! We also find 
that the cosmology cut restricts m t < 130(140) GeV for [i > 0(p < 0) in this range of 
nig. Therefore, the radiative correction contribution to m\ (oc m^) will be modest in this 
range of m g . This explains the depletion of points for rrih ^ 80 GeV in Fig. 5 and leads 



to the mass relationships in Eqs. fl5.4f) and ( |5 . 5| ) 



Recent studies of Higgs searches relevant to LEPII have shown that large regions of 
the parameter space which determine m^, including radiative corrections, can be explored 
|2^,^3|. These studies make simplifying assumptions regarding the many parameters at 



low energy, as well as choosing fixed values of mj. Nonetheless, we expect a certain level 
of quantitative agreement with these generic analyses. Important conclusions tend to be 
unanimous; if tan/3 ^ 5, then values of rrih up to rrih ~ 80 GeV can be explored (see 
e.g., Figs. 2a,8a in ||22|1 ). This is precisely the constraint on tan/3 that is realized in our 
analysis due to the stringent proton decay cuts. As we have discussed, the regions in 
parameter space where rrih ^ 80 GeV result in the constraint m + ^ 90 (110) GeV for 
«3 = 0.113 (0.120), and both signs of \i. Early model-dependent analyses of e + e~ Wlv\ 
indicate that ^ 100 GeV could be explored at LEPII [|4j. A more careful upper 



limit would require a detailed calculation, but these results are encouraging nonetheless. 
Therefore, if LEPII does not see the lightest Higgs, it has a good chance of seeing the 
lightest chargino instead, or viceversa. 



2 Note that Fig. 5 cannot be used to establish a lower bound on rrih since values of 1 < tan f3 < 2 
have not been considered. This is unlike the case for the quoted upper bound in Eq. ( |5.3| ). 
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6. Conclusions 

The most direct and pervasive evidence for unified models would be the observation 
of nucleon decay. In fact, this kind of test of a unified model has the very appealing 
property of involving the physics of both low-mass and high-mass particles in the theory, 
and as such should be able to discriminate among the various competing unified models 
at hand. In the specific case of the minimal SU(5) supergravity model, we have shown 
that under sensible assumptions, the range of proton lifetimes still to be probed by the 
next round of proton decay experiments is finite and completely accessible. Besides the 
various uncertainties on the measured 'constants' which enter the proton lifetime formula, 
two assumptions are key to the results presented in this paper: (i) the upper bound on 
the Higgs triplet mass Mh < SMjj, and (ii) the assumed upper bound on the squark 
and gluino masses of 1 TeV. Relaxing any of these assumptions can suppress the proton 
lifetime to acceptable values. However, we believe that if these quite sensible assumptions 
had to be relaxed to ensure compatibility with experimental requirements, then much of 
the motivation to consider this model would fade away and in effect the model would 
become 'sociologically' excluded. Another important variable in this study is the value of 
«3, larger values of which tend to open up the parameter space. We have considered the 
case of «3 = 0.113 and 0.120. A better defined value for this quantity will go a long way 
in nailing the specific predictions of the model. 

The recent precise observation by the COBE satellite of minute anisotropies in the 
cosmic microwave background radiation [^S[ is a dramatic reminder that precision cosmol- 
ogy has come to age, and that a unified model should be judged seriously by its consistency 
with cosmological constraints. In the present model we have shown that a suitable value 
for the cosmological relic density of the lightest neutralino is another powerful constraint 
on the model parameters. In effect it reduces the size of the allowed parameter space by a 
factor of 6. 

The proton decay cut singles out a small region of the parameter space where a 
number of fairly accurate mass correlations and mass bounds exists. We find for the one- 
loop corrected lightest Higgs boson mass rrih ^ 110 GeV, independently of the sign of 
the value of as, or the cosmological constraint. Since ^> Mw holds for most of the 
allowed parameter space, it also follows that m x ~ \ m x% ~ O-l^m^ and m x o « m x +- 
The imposed upper bound rriq < 1 TeV cuts off the value of the gluino mass for a given 
£o value and gives upper bounds on the light particle masses: m x < 85 (115) GeV and 
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to n + < 165 (225) GeV for a 3 = 0.113 (0.120). The effect of the cosmological cut on these 
predictions is negligible. However, this cut does remove a large portion of the allowed points 
and has a dramatic effect on the correlation between to/j, and m x + (or m x o , m x since they 
are all related), such that at least one of these particles is quite likely to be observable at 
LEPII. To recapitulate, of the 184,470 points in parameter space which we examined, only 
~ 1/240 satisfy the standard constraints, the proton decay bound, and the cosmological 
requirement, leaving only a rather restricted set of points to be put to experimental test. 
In sum, it should not be long before we could start discriminating among the various 
supersymmetric unified models. 
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Figure Captions 

Fig. 1. Scatter plot of the proton lifetime r p = r(p — >• u^ jT K + ) versus the gluino mass for 
the hypercube of the parameter space explored. The current experimental lower 
bound is Tp Xp = 1 x 10 32 y. The various 'branches' correspond to fixed values 
of £o as indicated (the labelling applies to all four windows). The bottom row 
includes the cosmo logical constraint. The upper bound on m g follows from the 
requirement rriq < 1 TeV. 

Fig. 2. The calculated neutralino relic density O x /iq versus the neutralino mass. Note 
the effect of the s-channel Z-pole at m x ~ \Mz- Only about 1/6 of the points 
have fi x /io < 1- 

Fig. 3. The calculated value of \i versus m g . The cosmological constraint is enforced in 
the bottom row. Above the solid lines (labelled B) the neutralino is a nearly 
pure bino eigenstate, above the line \n\ = M 2 = 0.3m g it is mostly a gaugino 
eigenstate. The small regions in the bottom row plots which extend up into the 
pure bino regions correspond to m x masses near \mz and \rrih- 

Fig. 4. The lightest neutralino mass (m x ) versus the gluino mass (top row) and the 
second-to-lightest neutralino mass (m x o ) versus the lightest chargino mass ( m x + ) 
(bottom row). These plots show the anticipated approximate mass correlations 
due to the proton decay constraint, m x ~ 0.15m g , m x o pa ~ 2m x . The 

cosmological constraint (not enforced in this figure) simply depletes the point 
density, without affecting the range of the particle masses. 

Fig. 5. The one-loop corrected lightest Higgs boson mass (m/J versus the lightest 
chargino mass. If the cosmological constraint is enforced (bottom row), then 
for nih ^ 80 GeV the lightest chargino is below 110 GeV. Thus, at least one of 
these particles is likely to be discovered at LEPII. A significant number of the 
allowed points correspond to m x masses near \Mz and \rrih. 
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